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Abstract We present an analysis of the diffuse soft X-ray emission from the nuclear region 
of M51 combining both XMM-Newton RGS and Chandra data. Most of the RGS spectrum 
of M51 can be fitted with a thermal model with a temperature of ^ 0.5 keV except for the 
O VII triplet, which is forbidden-line dominated. The Fe L-shell lines peak around the south¬ 
ern cloud, where the O viii and N vii Lya lines also peak. In contrast, the peak of the O vii 
forbidden line is about 10" offset from that of the other lines, indicating that it is from a spa¬ 
tially distinct component. The spatial distribution of the O vii triplet mapped by the Chandra 
data shows that most of the O vii triplet flux is located at faint regions near edges, instead of 
the southern cloud where other lines peak. This distribution of the O vii triplet is inconsistent 
with the photoionization model. Other mechanisms that could produce the anomalous O vii 
triplet, including a recombining plasma and charge exchange X-ray emission, are discussed. 

Key words: atomic processes - plasmas - ISM: jets and outflows - galaxies: Seyfert - 
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1 INTRODUCTION 


The Whirlpool galax y M51 (NGC 5194) is classifie d as a LINER or Seyfert 2 galaxy from studies of optical 
emission lines (e.g.. lStaufferlll982 : Ho et al. 1997 ). The close distance of M51 (7.8 Mpc, the mean value 


taken from NEE0 database) allows detailed studies of its nuclear activity. Radio observations show a bipolar 
outfl ow, comprising a southern cloud and a northe rn loop, which is also seen in the optical emission line 
map ( Eord et al.ll985 : Crane & van der Hulslll992 ). The optical outflow with velocities as high as 1500 km 
s“^ has been reported by Cecil ( 1988h . who also found that the arcuate radio emission of the southern cloud 


lies inside the bright optical emission line region, indicating a bow shock caused by the radio jet emanating 
from the nucleus. 

The interaction between the radio jet/outflow and the surrounding gas represents an important feed- 



X-ray morphology of M51 is similar to that seen in radio and optical observations. The X-ray spectra 
of the southern cloud and northern loop are similar and can be fitted by a thermal model of temperature 
~ 0.55 keV. The X-ray morphology and spectra suggest that the nuclear X-ray emitting gas of M51 are 
shock-heated by the bipolar radio outflow from the nucleus. Both the X-ray poi nt sources and the extended 
X-ray emission f rom the disk of M51 haye be en studied nsin^ Chandra (e.g., Terashima & WilsonI 2004 ; 


Tvler et aTll2004 1 and XMM-Newton data (e.g., Dewangan et aDl2005 : Owen & Warwickll2009 1. 
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Table 1 List of XMM-Newton RGS observations of M51 


ObsID 

Itot (ks) 

tea (ks) 

Obs time 

P.A. (degrees) 

0212480801 

49 

26 

2005-07-01 

294 

0303420101 

54 

35 

2006-05-20 

326 

0303420201 

37 

25 

2006-05-24 

323 

0677980701 

13 

13 

2011-06-07 

312 


Note: ttot is the total exposure, is the useful exposure after removing periods of flares, and P.A. is the position angle. 


All these X-ray studies are based on CCD data, for which the spectral resolution is arou nd 10 at 1 keV. In 


contr ast, the Reflection Grating Spectrometers (RGS) aboard the XMM-Newton telescope (Iden Herder et al. 


200 Ih have a better spectral resolution for moderately extended sources due to their large dispersion power. 
The RGS spectral resolution is ^ 0.140 A, where 9 (in units of arcmin) is the angular extent of the source 
along the dispersion direction. For M51, the spatial extent of the nuclear X-ray emitting region is about 
0.5', which corresponds to a spectral resolving power X/5x 150 at 15 A given the spatial resolution 
of XMM-Newton. With this resolution many emission lines, especially the He-like O vii triplet which is a 
diagnostic tool for the thermal state of plasma, are well resolved. 

The O VII triplet consists of a resonance line, two inter-combination lines and a forbidden line. For an 
optically-thin thermal plasma in ionizat i on eq uilibrium, the electron collisional excitation is efficient and 

(l2012h . we studied nine nearby star-forming galaxies (including 


favors the resonance line. In Liu et al 


M51) and found that the forbidden lines of their O vii triplets are comparable to or even stronger than the 
resonance lines. We proposed the charge-exchange process between highly ionized ions and neutral species 
as a possible explanation. The charge-exchange captured elect rons of the recip ient ions are in excited states 
and their downward cascading favors the forbidden line (e.g. Dennerll 2O10li . In the special case of M51, 
however, the explanation is complicated by the presence of the central low-luminosity AGN and the radio 
jet-driven outflow. 

In lLiu et al.l (120121) . only the line ratio of the O vii triplet of M51 was studied. In this paper we provide 
a more detailed study of the XMM-Newton RGS spectrum of M51. In particular, we study the spatial dis¬ 
tribution of emission lines along the cross-dispersion direction of RGS. As shown in Sections, the spatial 
distribution of the O vii triplet is different from other lines. Given this situation, the Chandra image of M51 
with sub-arcsec angular resolution is helpful in revealing the O vii triplet distribution. Thus we also analyze 
the archival Chandra data of M51 with an accumulated exposur e time of 700 ks, which is m uch deeper than 
the short exposure data (^ 15 ks) used in previous studies (e.g. Terashima & Wils^l2001 ). 

We describe the observational data in Section 2 and present the results of XMM-Newton and Chandra 
data in Section 3 and Section 4, respectively. Discussion of the results is covered in Section 5. Throughout 
the paper, the errors quoted are for the 90% confidence level. At a distance of 7.8 Mpc, 1 arcsec corresponds 
to 38 parsec. 


2 OBSERVATIONAL DATA 

The two RGSs on-board the XMM-Newton telescope are slit-less dispersive spectrometers, and photons 
from extended sources are recorded on CCD detectors with the dispersion angle and ID spatial information 
along the cross-dispersion direction. Because the dispersion directions are different for different observa¬ 
tions, to ensure that the dispersed spectra are from similar spatial regions, we use four archival datasets of 
XMM-Newton RGS observations of M51 as listed in Table 1. The maximum difference of position angles 
between them is about 30 degrees. The total effective exposure time is ^ 100 ks after removing intense flare 
periods. 

The most recent version of the Science Analysis System (SAS 14.0) of XMM-Newton is used for the 
reduction of photon events. As stated above, the RGS spectrum is broadened by the spatial extent of the 
source along the dispersion direction. To produce a broadened redistribution matrix file (RMF), we con¬ 
volve the RMF produced by the SAS tool rgsproc with a Chandra image of M51 (0.4-2 keV) using the 
rgsrmfsmooth tool written by Andrew Rasmussen. The broadened RMF is calculated separately for each 
dataset. For background subtraction, we use the model background generated by rgsproc based on the flux 
of CCD 9 beyond 1' from the on-axis position. Because M51 lies at ~ 1' away from the on-axis position 
for most observations listed in Table 1, the model background is over-predicted. Indeed, the model back¬ 
ground exceeds the observed spectrum at A > 28 A. Thus we scale down the model background by a factor 
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Dispersion (A[A]) 

Fig.l Cross-dispersion vs dispersion CCD image of M51 combining both RGS1 and RGS2 
data of all four observations. The emission lines are seen as vertical contours and are well re¬ 
solved. The color bar is for the total event counts on linear scale. The zero point along the cross¬ 
dispersion direction of ObsID 0303420201 is adopted. 



Fig. 2 Chandra counts image of M51 within 0.4 — 2 keV (logarithmically scaled). The overlaid 
ra dio contours (1.5, 3,6 and 12 in units of 10 pJy) are from VLA observations at 6 cm wavelength 
bv ICrane & van der Hulslj (Il992h . The two solid lines separated by 30" indicate the RGS disper¬ 
sion direction of ObsID 0303420201, while the arrow indicates the positive cross-dispersion 
direction. The white cross marks the position of the nucleus of M51, which corresponds to a 
cross-dispersion distance of 52" for ObsID 0303420201. 


of 0.6 based on the event distribution of CCD9. Since we focus on emission lines, the uncertainty in the 
background subtraction will not affect our major results. 

Because CCD 4 of RGS2 covering 20 — 24 A failed early in the mission, the data of the O vii triplet 
around 22 A are only from RGSl. Similarly, CCD 7 of RGSl also failed and the data within 10.6 — 13.8 
A of RGS 1 are missing. Figure 1 presents the CCD image of the XMM-Newton RGS data of M51 combining 
both RGS 1 and RGS2 for all four observations. The offsets between RGS 1 and RGS2 and between different 
observations have been corrected. The figure shows clearly the well-resolved emission lines. 

We use seven Chandra datasets with ObsID numbers of 13812, 13813, 13814, 13815, 13816, 15496, 
15553, as observed by PI Kip Kuntz. After removing the flare periods, the total effective exposure time 
is about 700 ks. The datasets are analyzed with CIAO (version 4.6) following standard procedures. As an 
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Fig. 3 Spatial distribution of the emission lines of Fe xvii at 15 A, O viii Lya at 19 A, and O vii 
forbidden line at 22.1 A. For clarity, the profile of the Ovii forbidden line at 22.1 A has been 
multiplied by a factor of 4. 



Wavelength (A) 


Fig. 4 Background-subtracted XMM-Newton RGS spectrum of M51 combining both RGS1 and 
RGS2 data of all four observations. Red histograms are the fitted thermal model with four 
Gaussians. x is the difference between data and model divided by the error. 


illustration. Figure 2 presents the counts image of M51 within 0.4-2 ke V created by merging all the seven 
datasets. The over-plotted contours are from the radio 6 cm observation ( Crane & van der Hulstlll992 ). and 
show a close correspondence with the X-ray image. The XMM-Newton RGS dispersion direction is also 
illustrated in Figure 2. 


3 XMM-NEWTON RGS RESULTS 

3.1 Spatial distribution of X-ray emission lines 

First we study the spatial distribution of the well-resolved X-ray emission lines along the cross-dispersion 
direction of XMM-Newton RGS. The Fexvii lines at 15 and 17 A, the O viii and N vii Lya lines at 19 and 
24.8 A, and the O vii forbidden line at 22.1 A are studied. The wavelength regions adopted to calculate the 
photon counts are 14.8 - 15.5 A, 16.6 - 17.4 A, 18.6 - 19.4 A, 24.4 - 25.2 A, and 21.8 - 22.4 A, respectively. 
We find that the profiles of Fe xvii lines at 15 and 17 A are similar, as are the O viii and N vii Lya lines. For 
clarity, only the profiles of Fe xvii line at 15 A, O viii Lya line at 19 A, and O vii line at 22.1 A are plotted 
in Figure 3. We see that the bright Fexvii line at 15 A is centrally peaked around 50", with a FWHM (full 
width at half maximum) of ~ 20". This FWHM is similar to the spatial resolution of XMM-Newton RGS. 
Since the cross-dispersion position of 50" corresponds to the location of the southern cloud, it means that 
the emission lines of Fe xvii are mainly from the compact southern cloud. 

The O VIII Lya line at 19 A is also peaked around 50", the same as the Fe xvii lines. The spatial distribu¬ 
tion of the Ovii forbidden line at 22.1 A, however, is different from those of the other lines. Its maximum 
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Table 2 Fitting results 


T (keV) 

N 

O 

Ne 

Fe 

O yiir 0 yii / 

Nyir 

N Ylf xi 

0.52 ±0.02 

4.4 ± 1.1 

0.6 ±0.1 

0.6 ±0.2 

0.3 ±0.1 

1.3 ±0.7 4.1 ± 1.0 

2.2 ±0.9 

2.5 ± 1.2 0.91 


Note: the abundances are relative to the solar values; OVII,.,/ and NVI^ y (in units of 10 ® photons s ^ cm ai'c the line 
intensities of O VII resonance line, O VII forbidden line, N VI resonance line, and N VI forbidden line, respectively. 


is around 60", 10" offset from the peaks of other lines. This indicates that the O vii forbidden line is from a 
spatially distinct component, compared with the other lines. 

3.2 X-ray spectrum of M51 


Since the X-ray emission lines have different spatial distributions, it is optimal to extract spectra within 
different spatial regions. However, the difference of the scale of ^ 10" is too small to allow such spectral 
extractions from XMM-Newton RGS. Thus we only extract one spectrum centered on the nucleus of M51 
with a cross-dispersion width about 1'. The background-subtracted spectrum of M51 combined from both 
RGS 1 and RGS2 data of all four observations is plotted in Figure 4. 

We see that the spectrum of M51 is dominated by emission lines. The O vii triplet is dominated by the 
forbidden line at 22.1 A and the N vi triplet also has a forbidden line comparable with its resonance line. The 
bright Fe L-shell lines around 15 - 17 A indicate that they are not from a photoionized plasm a, for which the 
Fe L-shell lines are expected to be much weaker (e.g.. Kallman et al.ll996 : Sako et al .l2nnnh. They are likel y 
due to a collisionally ionized plasma as suggested in preyious studies (e.g., Terashima & Wilson 200 ih. 
Thus we fit an optically thin collisional-ionization-equilibrium (CIE) thermal model lyapec. iFoster et al. 


2012h to the obseryed RGS spectrum of M51. Only the abundances of the elements of Ne, Fe, O, and N, 
which show bright emission lin es in the spect rum of M51, are allowed to yary. The abundances of other 
elements are set to solar yalues ( Lodderj2003 ). 


On the other hand, the forbidden-line dominated O vii triplet of M51 shows a spatial distribution differ¬ 
ent from other lines. As stated in Section 1, it is impossible to explain such an anomalous triplet by a thermal 
CIE model. Thus we add two Gaussians to represent the Ovii resonance and forbidden lines. The inter¬ 
combination line is too weak to be fitted. Similarly, we add another two Gaussian s to represent the N yi res¬ 
onanc e and forbidden lines. The model is subject to an absorption model of wabs ( Morrison & McCammonI 
1983 1 with a foreground neutral hydrogen column density of 2 x 10^°cm“^ ( Kalberla et alJ ^OSh . The 
fitted results are listed in Table 2 and oyer-plotted in Eigure 4. 

The model proyides a reason able fit to the obsery e d spec trum of M51. The fitted temperature is 0.52 
keV, similar to that obtained by iTerashima & WilsonI (1200111 . Howeyer, the fitted abundances are higher 
than their yalues, which are only around 0.1 solar yalue. The low abundances they obtained are most likely 
due to the poor spectral resolution of Chandra data. We note that the fitted O and N abundances are not 
affected by the Gaussian fitting to the corresponding O vii and N vi triplets, since they are most likely due to 
a mechanism different from CIE models and at the fitted temperature, the emissiyity of O yii and N yi triplets 
of the CIE model is weak. The residuals around the Ee and O lines are likely due to the simple conyolution 
modeling of the RME There are residuals around 13.5 A, indicating that the Ne x triplet is under-estimated. 
The residuals around 18.6 A are likely due to the O yii He/3 line. 

The temperature dia gnostic G ratio of the He-like triplet for a thermal CIE plasma is defined as 
( Gabriel & Jordanlll969h 


G = 


f+i 


( 1 ) 


where f, i, and r represents the intensity of the forbidden, inter-combination, and resonance li nes, respec 


tiyely. If we assume the intensity of the inter-combination line is 1/4 of that of the forbidden line ((Smith et al. 


200 Ih . the fitted G ratio of the O vii triplet is 3.9 ± 2.3. Similarly, the fitted G ratio of the N vi triplet is 
1.4 ± 0.8. The G ratio of Ovii triplet shows that the Oyii triplet is not due to CIE thermal plasmas, for 
which the G ratios are generally smaller than 1. The G ratio of N vi triplet is barely consistent with a CIE 
model around 10® K. 
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Fig. 5 Left: Fluxes of the O vii triplet obtained with the Chandra data of M51 using the contour 
binning method. Right: The binned counts map of M51 within 0.4 - 2 keV. The contours of levels 
of 20, 50, 80, 200, and 550 for unbinned map are plotted in both panels. The white cross marks 
the nucleus of M51. 


4 Chandra RESULTS 


As shown in Sections. 1, the peak of the spatial profile of the O vii forbidden line is offset from those of the 
other lines by ~ 10". This scale can be easily resolved with Chandra data. Thus we analyze the Chandra 
data of M51 in this section. _ 

We use the contour binning method ( Sandersll2006 ) to divide the nuclear X-ray emitting region of M51 
into 34 different bins with a minimum signal to noise ratio of 40. The X-ray spectrum is extracted from 
each bin with a background spectrum taken from a source-free region outside the nuclear region. We fit a 
thermal CIE model with variable O, Ne, Mg, and Fe abundances and limit the fitting range to 0.4 — 1.5 
keV. We include a Gaussian (with a center of 21.85 A and a width of 0.25 A) to represent the O vii triplet. 
Adding a Gaussian will improve the by 4.3 on average. The improvement of is 8.1 for the spectra 
with the Gaussian line fluxes larger than 5 x 10”^ photons s“^ cm“^. We also tested the fitting method by 
fitting mock spectra with a line flux of 1 x 10“® photons s“^ cm“^, and found that about 75% of the fitted 
fluxes are within 3 sigma of the true value. 

Due to the limited spectral resolution of Chandra data, we do not expect the fitted O vii triplet fluxes to 
be as accurate as the measurement by the XMM-Newton RGS data. Nevertheless, it is likely to provide a 
qualitative distribution of the O vii triplet, especially for bins with high O vii triplet fluxes. The fitted fluxes 
of the Ovii triplet from the Chandra data are plotted in the left panel of Figure 5. For comparison, the 
binned counts image is plotted in the right panel. We see that the fluxes of the O vii triplet are generally 
higher for faint bins near the edges. It shows that the fitted Ovii triplet has a more extended distribution 
than that of the total emission within 0.4 — 2 keV, which is centered on the southern cloud. The northern 
regions have more O vii triplet flux than the southern part. This is consistent with the XMM-Newton result. 
Another interesting feature to note is that the southern cloud seems to be enclosed by an arc of O vii triplet 
in the outward direction. 


5 DISCUSSION 

We have analyzed the XMM-Newton RGS spectrum of M51, the emission lines of which are well resolved. 
Most of the spectrum of M51 can be fitted with a thermal CIE model of temperature ~ 0.5 keV, except for 
the Ovii triplet, which is forbidden-line dominated. The Eexvii lines at 15 and 17 A are centrally peaked 
around the cross-dispersion position of 50" with a EWHM of 20". This indicates that they are mainly from 
the compact southern cloud. The Oviii Lya line at 19 A and Nvii Lya line at 25 A also peak around the 
same position as the Eexvii lines. In contrast, the peak of the O vii forbidden line is about 10" offset from 
the other lines. The Ovii triplet map obtained with the Chandra data shows that most of the fluxes of the 
O VII triplet are located at faint regions, instead of the southern cloud where the other lines peak around. 
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As stated in Section 1, the forbidden-line-dominated O vii triplet of M51 cannot be due to a thermal 
CIE plasma, and we have proposed the charge-exchange process as a possible explanation. Nevertheless, 
the existence of low-luminosity AGN and the radio jet-driven outflow in the nuclear region of M51 makes 
other explanations possible, including photoionization and non-equilibrium-ionization plasmas. Below we 
discuss them in turn. 

For a photoionized plasma, the emission is dominated by recombination, which also favors the forbidden 
line. One characteristic feature of the photoionization model is the spatial profile. As the central AGN is 
the ionizing source, the ionizing flux will decline as r~^ and photoioni zation is most impor tant for regions 
close to AGN. This has been illustrated by the optical study of M51 (IBradlev et al.ll2004ll . which shows 
that photoionization is dominant within the inner region (r < 1") and the shock model is preferred outside 
(r ~ 2.5"). The mapped morphology of the Ovii triplet from Chandra data shows no fluxes around the 
nucleus and is most prominent at the outer regions. This is inconsistent with the photoionization model. 

All the optical, radio and X-ray studies of M51 support a scenario in which the nuclear gas of 


M51 is shock-heated by the radio outflow emanating from the nucleus (ICecillIl988t iFord et al.l 11985 


Terashima & Wils^l2001 ). In this situation, the non-equilibrium-ionization (NFI) plasma is another possi 


ble explanation for the anomalous O vii triplet. If a plasma is shock-heated suddenly, the ionization process 
lags behind heating, and will result in an ionizing plasma. The inner-shell coll isional ionizat ion of 
ions can lead to excited 0®+ ions and enhance the forbidden line emission (e.g., Fiedahll 1999h . However, 
this on ly happens in a short t ransition period when the 0®“*" fraction is non-negligible. Using the Sedov 
model ( Borkowski et al.ll2001 ), we And that the G ratio of the O vii triplet is around 3 at an ionization time 
net ^ 3.5 X 10® cm“®s. This timescale is too short compared with the spatial extent we studied. A more 
likely explanation is a recombining NFI plasma. 

When a blast wave expands into a rarefied medium, a recombining plasma could be produced due to 
rapid adiabatic cooling, as prop osed to explain the radi ative recombination continuum features in some su¬ 
pernova remnants (SNRs) (e.g., Yamaguchi et al.ll2012 ). The spatial distribution of the O vii triplet mapped 
by the Chandra data is consistent with this scenario. A unique test of a recombining plasma is the associ¬ 
ated radiative recombination continuum. However, due to the presence of strong Fe lines, it is hard to tell 
whether there is a radiative recombination continuum of O vii (at 16.8 A) in the spectrum of M51. 

On the other hand, to explain the anomalous O vii triplet with the charge-exchange process, the interact 


ing in terfaces b etween the highly ionized gas and neutral species are needed. The obseryed HI (IWalter et al 


2008h and CO (IKoda et al.N201 Ih maps of M51 are well correla ted with the spirals. CO is also detected 


on 1-2 arcsec scales around the nucleus (IMatsushita et al.ll2007h . Currently there is no much information 
of neutral materials on 10" scales of M51. Further obseryations are needed to test the existence of neutral 
materials on 10" scale. Different from the recombining plasma, the charge-exchange model shows no ra- 
diatiye recomb ination continuum, but ha s enhanced emission lines from preferred energy leyels (n = 4 — 7 
for OVII) (e.g. iBeiersdorfer et al.ll2003h . Future spatially-resolyed high-resolution X-ray spectroscopy by 
the soft X-ray spectrometer of Astro-H will be able to test the existence of the radiatiye recombination 
continuum of M51. 
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